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Abstract. Spitzer Space Telescope observations of a point-like source, L1014-IRS, close to the dust peak of the low-mass 
dense core L1014 have questioned its starless nature. The presence of an object with colors expected for an embedded protostar 
makes L1014-IRS the lowest luminosity isolated protostar known, and an ideal target with which to test star formation theories 
at the low mass end. 

In order to study its molecular content and to search for the presence of a molecular outflow, we mapped L 10 14 in at least one 
transition of l2 CO, N 2 H + , HCO+, CS and of their isotopologues 13 CO, C 18 0, C 17 0, N 2 D + and H 13 CO + , using the Five College 
Radio Astronomy Observatory (FCRAO), the IRAM 30 meter antenna and the Caltech Submillimeter Observatory (CSO). 
The data show physical and chemical properties in L1014 typical of the less evolved starless cores: i.e. H 2 central density of 
a few 10 s molecules cm 4 , estimated mass of ~2 M Q , CO integrated depletion factor less than 10, _/V(N 2 H + )=: 6 x 10' 2 cm -2 , 
iV(N 2 D + )//V(N 2 H + ) equal to 10% and relatively broad N 2 H + (l-0) lines (0.35 km s _1 ). Infall signatures and significant velocity 
shifts between optically thick and optically thin tracers are not observed in the line profiles. 

No classical signatures of molecular outflow are found in the 12 CO and 13 CO observations. In particular, no high velocity wings 
are found, and no well-defined blue-red lobes of 12 CO emission are seen in the channel maps. If sensitive, higher resolution 
observations confirm the absence of an outflow on a smaller scale than probed by our observations, L1014-IRS would be the 
only protostellar object known to be formed without driving an outflow. 

Key words. ISM: clouds - ISM: evolution - ISM: individual(L1014) - ISM: molecules - ISM: jets and outflows - Stars: 
formation 



1. Introduction 

Starless cores are cold (~ 10-20 K) and dense (> 10 4 cnT 3 ) 
condensations of gas and dust in which no sign of the pres- 
ence of a central protostellar object has been found. Previously, 
the easiest way to determine if dust emission was associated 
with a protostellar object was to search for a source emit- 
ting at mid-infrared w avelengths in the IRAS catalogue (see, 
Beichma n et all [l986). This technique was obviously limited 
by the IRAS sensitivity. 

A clear example is rep r esent ed by IRAM 04191: although 
Andre, Motte & Bacmann ll 19991) found a clear sign of the 
presence of a Class object represented by a clear collimated 
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bipolar outflow departing from the core peak, IRAS does not 
show a point source towards the dust emission peak. Moreover, 
the CS line they observed towards the nucleus confirmed the 
presence of star formation activity showing a clear double 
peaked profile with the blue peak brighter than the red one. 
This spectral asymmetry ha s been recognized as an indicator of 
systematic i nward motion dSnell & Lorenl ll977tlZhoulll992t 
Tafalla et ail ll998tlLee. Mvers & Tafallalll999h . 



In the same fashion, we present here the c ase of L 1014. 
This core, listed as an opacity class 6 object in the lLvndsl(ll962l) 
catalogue, lies -0°25 below the galactic plane, has a line of 
sight velocity of ~4 km s _1 with respect to the local standard 
of rest and is projected just 10' south of another dark glob- 
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ule B362 1 . These cores can be seen in the Digital Sky Survey 
optical image presented in Figure ^ where we have overlaid 
visual extinction contours derived from near-infrared color ex- 
cesses of background st ars listed in the 2MASS catalo g using 
the NICE technique (e.g. Lada et allll994lAIves et allll998h . 
L1014 was included in a survey for infall asymme- 
try in starle ss cores performed observ ing CS(2-1) and 
N 2 H + (l-0) fc.ee. Mvers & Tafallai Il999h : given the weak 
detection in CS, no N 2 H + observations were attempted. 
The dust continuum from L1014 was first detected at 
260 u m using the NASA Kuiper Airborne Observatory 



rKAOtlKeenel7 l981). Subsequently a 850 /im emission map 
llVisser. Richer & Chandler! EoOl) taken with SCUBA showed 
a dust emission peak at (21 h 24 m 07.6, 49°59'02", J2000) coin- 
cident with the visual extinction peak. These authors then per- 
formed an unsuccessful search for high velocity gas, observing 
12 CO(2-l) at five points around the dust peak with a sensitivity 
of 0.3 K km s _1 in a 0.2 km s -1 channel. 

The Spitzer Space Telescope (hereafter Spitzer) observed 
L1014 in December 2003 as part of the Leg acy program 
"From Molecular Cores to Planet Forming Disks" dEvans et all 
120031) . Surprisingly, a strong (81.8+16 mJy) point source was 
dete cted with MIPS at 24 pun towards the center of L1014 
.Young et all l2004. being coinc i dent w ith both the 850 /im 
peak (|Vis ser. Rich er & Chandler! 120011) and 1 .2 mm peak 
.Young et alJl2 004: Kauffmann et al., in preparation). This ob- 
ject, referred to as L1014-IRS, is also clearly detected in each 
of the 4 IRAC bands (3.6, 4.5, 5.8 and 8.0 /zm) and at 70 /im 
with MIPS and, most importantly, shows a Spectral Energy 
Distribution (SED) compatible with an embedded protostar 
having an effective temperature of 700 + 300 K according to 
the IRAC data at A < 8 /im. Considering that L1014-IRS is 
the only source visible at 70 /mi in the 5'x5' field observed 
with the Spitzer and L1014 is the on ly starless core present in 
the same region, lYoung et all .2004 concluded that a chance 
alignment of the two sources was unlikely although not neg- 
ligible. Naively, dividing the Spitzer field of view in squared 
cells with 30" of side, the chance probability of having both 
L1014 and L1014-IRS in the same cell would be 1%. In this 
paper, although we will consider the chance alignment as an 
alternative explanation, we will devote more discussion to the 
hypothesis of association since the weight of the evidence 
supports this hypoth esis. A ssuming a distance to the core of 
200 pc, Youn g et all d2004 determine a luminosity of 0.09 L 
for L1014-IRS. From the bolometric temperature (Tb i= 50 K) 
and the ratio betwee n bolometric and s ub-millimeter luminos- 
ity (L ho i/L smm = 201. lYoung eTall ll2004 classified L1014-IRS 
as a Class protostar (see lAndra.ll994t for a review on Class 
properties). 

Outflows are an ubiquitous tra cer of protostellar activity 
llAndrel 1 19941 iRicher et all l2000h . and so, in an attempt to 
support the classification of L1014-IRS as a Class protostar, 
we deepened the search for high velocity 12 CO with a better 
combination of resolution and sensitivity than in past searches. 



We report here on observations of 12 CO, 13 CO, C ls O, C 17 0, 
CS, N 2 H + , N 2 D + , HCO + and H 13 CO + aimed at surveying the 
molecular content of L1014 and at studying its kinematical 
properties. 



2. Observations 

Observations were performed during March and April 2004, 
using the 32-pixel SEQUOIA focal plane array mounted on the 
Five College Radio Astronomy Observatory (FCRAO) 14 me- 
ter telescope. 

We mapped L1014 in CO(1-0), I3 CO(1-0), C I8 O(l-0), 
N 2 H + (l-0) and CS(2-1) using the on-the-fly position switch- 
ing mode. Adopted frequencies, telescope half power beam 
widths, system temperatures, channel spacings and size of the 
mapped area are in Table [2 Temperatures were converted in 
the main beam brightness scale according to the efficiencies 
tabulated in the FCRAO web page 2 . Data were reduced and 
convolved onto a 20" spa ced grid (circa Nyquist sampling) us - 
ing the otftool software (Hever. Narayanan & Brewer , 200 1) 3 
available at FCRAO and then analysed using CLASS, the 
spectral line analysis software of IRAM and Observatoire de 
Grenoble (see Bui sson et alll2002l for a manual). 
Angular offsets throughout this paper are relative to the 1.2- 
mm continuum peak position (21 h 24 m 07 s .5, 49°59'05", J2000), 
which is 4" south of L1014-IRS position. 

Complementary higher-resolution observations of 
N 2 H + (l-0), (3-2), N 2 D + (l-0), (2-1), (3-2) and C 17 0(1- 
0) were obtained in August 2004 using the IRAM-30m 
telescope. 

The observations consisted in a 5-pointings cross spaced 
by 20" around the 1.2-mm peak position for N 2 H + (l-0), 
N 2 D + (l-0) and N 2 D + (3-2), while N 2 H + (3-2) and N 2 D + (2-l) 
were observed only at the (0,0), (20,0) and (-20,0) and 
C 17 O(l-0) only at the peak position. These data were taken 
in frequency switching mode using the facility autocorrelator 
as the backend. Information about the telescope parameters 
(channel spacing, HPBW, system temperatures) are reported 
in Table |2 Intensities were converted in the main beam 
brightness scale using the efficiencies reported in the IRAM 
web site 4 . 

Finally, we used the Caltech Submillimeter Observatory in 
August 2004 to obtain HCO + (3-2), H I3 CO + (3-2) and HCN(3- 
2) spectra towards the central position and a small map in 
C ls O(2-l) and C 17 0(2-l). Spectra were obtained in position 
switching mode. Although the single pointing observations 
were performed 6" south of L1014-IRS, the ~30" beam al- 
lows us to consider them as relative to the infrared source posi- 
tion. Conversion factor to the main beam temperature (t] 210GHj , = 
0.64) was evaluated from observations of planets. Other infor- 
mations for the telescope set-up are reported in Table[2] 

The pointing accuracy for all our observations was mea- 
sured to be as 4 - 6". 



1 We note that B362 and L1014 were called "L1014-1" and "L1014- 
L in Lee & Mversl ll999h. In the present paper we keep the o riginal 
names as in Barnard. Frost & Calvertl J 19271) and lLvnds|]l962T) . 



http://www-astro.phast.umass.edU/~fcrao/observer/statusl4m.html#ANTENNA 

http://www-astro.phast.umass.edu/~fcrao/library/manuals/otfmanual.html 

http://www.iram.fr/IRAMES/telescope/telescopeSummary/telescope_summary.htr 
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3. Results 

3.1. Maps 

In Figure |2] we present integrated intensity emission maps 
of L1014 (Vlsr =* 4.2 km s" 1 ) in N 2 H + (l-0), CS(2-1) 
and C 18 O(l-0) take n with FCRAO to gether with the 1.2- 
mm continuum from lYoung et al.l J2004I) . the 5-points map of 
N 2 H + (l-0) taken with the IRAM-30m and the 13-points map 
of C I8 0(2-1) from CSO. The half maximum contour of each 
map is rendered in white. The white cross shows the posi- 
tion of L1014-IRS, the candidate pr otostel lar object embedded 
in L1014 according to lYoung et al.l ( | 2004). The area observed 
with the FCRAO is 10'xlO' in size, but only the inner 8'x8' 
region where emission was detected is shown in Figure|2] 

N2H + (l-0) emission, which comes from the high density 
nucleus of the core, is very concentrated (FWHM ~70", in 
the FCRAO observations) and peaks ~10" west of the dust 
emission (within the uncertainties due to the 54" FCRAO 
beam width and the Nyquist sampling of the observations). 
Since the emitting region is comparable in size with the tele- 
scope beam some dilution is very likely. In fact, IRAM-30m 
observations of N2H + (l-0) with a 26" beam showed spectra 
3 times brighter than those observed by FCRAO towards 
the dust continuum peak, although still the emission does 
not seem to be fully resolved. In these higher resolution 
observations, the integrated intensity peak falls 5" south 
of L1014-IRS. We note that the equivalent radius of the 
N 2 H + (l-0) emission from FCRAO (=0.034 pc, evaluated as 
the square root of the area within the 50 % contour divided by 
n) is a mong the smallest radii found inlLee. Myers & Tafalla 
J200ll) and ICaselli. Benson. Mvers & Tafallal fc002t . which 
observed 63 starless cores altogether using FCRAO and 
evaluated the equivalent radius with the same tech- 
nique. The integrated intensity of N2H + (l-0) at peak 
observed with FCRAO is also low er (a factor of 2) 
than any other core reported in |Lee, M vers & Ta falla 
fcOO ill and ICaselli. Benson. Mvers & Tafallal d2002l) . 
although th e beam dilution should be k e pt in mind. 
Following Casel li. Benson. Mvers & Tafallal J2002I) . we 
calculated the virial mass of L1014 from M v i r (M ) = 
210 r (pc) Av 2 (km 2 s" 2 ), where r is the typical radius 
of the N7H + emission from FCRAO and Av 2 is the 
velocity dispersion of the m e an m ass molecule (see 
ICaselli. Benson. Mvers & Tafallal 120021 for the defini- 
tion). In this way we obtain ed a virial mass of 2.1 M , again, 
among the smallest in the ICaselli. Benson. Mvers & Ta falla 
(2002) sample. This estimate is in agreement with the mass 
evaluated from the extinction map in Figure ^ In fact, using a 
conversion factor of 1 X 10 21 molecules cm" 2 per magn itude 
of visual extinction (Bohlin, Savage, & Drake, 1978), we 
obtained a mass of 1.5 M above the Ay = 5 magnitudes 
contour. 

The CS(2-1) emission is much more extended than the 1.2- 
mm continuum and N2H + (l-0), as expected du e to its lower 
critica l density and possible depletion effects jTafalla et all 
120021) . The integrated intensity peak lies ~30" west of the dust 



emission, although the same caveats about the limited resolu- 
tion and sampling apply here as in the N2H + map. 

C 18 O(l-0) observations reveal more structure than previ- 
ous observations, showing a secondary peak at (-200", -200") 
from the reference position. Note that the 1.2-mm map did 
not cover that area. Also the C 18 O(l-0) integrated intensity 
peaks west of the continuum, although the peak is less con- 
strained than the other tracers. The flatness of the C ls O map 
and the offset of its peak with respect to the dust are typical 
features of starless cores (Tafalla et al., 2002), where CO (and 
its isotopologues) freeze-out onto dust grains in the high den- 
sity nucle us. To t est if t his is the situation in L1014, we fol- 
low ICaselli et alJ ll2002bl) evaluating the degree of CO deple- 
tion as the r atio of the canonical CO abundanc e ([CO]/[H 2 ]= 
9.5 x 10" 5 , Frerking. Langer& WilsorJ Il982l) and the ob- 
served abundance derived fro m the ratio of C 18 Q and H2 col- 
umn densitie s (for details see ICaselli etall E002b). Following 
ICrapsi et alJ (^005), the CO column density was derived un- 
der the constant excitation temperature approximation (CTEX, 
see eq. A4 in ICaselli etall Eo02b) assuming that C 18 emis- 
sion is optically thin. In this way we derived a C ls O col- 
umn density of 7.8 x 10 14 cm" 2 at the dust peak position 
which increases to 8.2 x 10 14 cm" 2 at the C 18 peak posi- 
tion (-20", 20"). The CO column density can then be inferred 
from the local inte rstellar medium relati ve abundance value 
(rCOl/rC 18 Ol=560: IWilson & Roodll994T) . H 2 column density 
was inferred from the 1 .2-mm emission assuming constant dust 
temperature (10 K) and emissivity (K] ? mm = 0.005 cm 2 g -1 ) 
similarly to IWard-Thompson. Motte & Andre) |l999), obtain- 
ing N(H2) — 4 x 10 22 cm" 2 at the dust peak position and 
N(H 2 ) = 2.2 x 10 22 cm" 2 at the C ls O peak position. We thus 
found a depletion factor around 9 at the dust peak, whereas it 
drops to 5 at the C ls O peak. These values are typical of starless 
cores with moderate chemical evolution (e.g. L 1495, L1498 , 
L492 ) with densities of about a few 10 5 cm" 3 llCrapsi et all 
2005). 

The 30" map of C 18 0(2-l) obtained at CSO is shown in an 
inset of Figure^] Note that the error on the integrated intensity 
for the points outside the peak is 2.5 times the value reported 
in Tabled Similarly to the N2H + high-resolution observations, 
the C 18 0(2-l) from CSO seems to peak south of L1014-IRS. 

While the 1.2-mm emission peak in L1014 (23 mJy (11" 
beam)" 1 ) is comparable in intensity with the maj o rity o f 
N2FF emitting starless cores (see e.g ICrapsi et all [2005), 
it is relatively weak if compared with high-density, evolved 
pre-stellar cores such as LI 544 (> 60 mJy (11" beam)" 1 , 
Ward -Thompson, Motte & Andre 1999). Consequently, the H2 
central volume density inferred from it is 6 times less (2.5 x 
10 5 cm" 3 in L1014 vs. 1.4 x 10 6 cm" 3 in L1544) when eval- 
uated with the sam e technique and physical parameters as in 
iTafalla et all J2002I) . This technique derives the density profile 
of a spherically symmetric core which best fit the observed 1 .2- 
mm continuum. In our calculations we adopted a density profile 
model of the form n(H2) = «o/(l + f/fo) a , finding, in addition 
to the central density, the radius of the "flat" region ~20" and 
the asymptotic power law index ~2 .7. Note that the central H2 
density value given in lYoungetalJ 12004) (1.5 X 10 5 cm 3 ) was 
evaluated using a dust opacity equal to 0.0104 cm 2 g" 1 . 
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3.2. Spectra towards dust peak 

The spectra observed towards the dust emission peak with 
FCRAO are shown in Figure [3] while those observed with 
IRAM-30m are shown in Figure |4] and those taken at CSO in 
Figure |5] Gaussian fits were performed using the CLASS fit- 
ting procedure and results are presented in Table[3] In the case 
of lines with hyperfine structure, all the hyperfine components 
were simultaneously fit. This kind of fit also furnishes the opac- 
ity of the line in the presence of high S/N spectra. The N2H + (1- 
0) line widths (0.35 km s _1 ) are slightly broader than the proto- 
type of evolved starless core L1544 (0.31 km s^. lCaselli et all 
2002a), suggesting a combination of nonthermal and system- 
atic motions in the inner nucleus of L1014. Using the CTEX 
approximation on the higher resolution (IRAM-30m) data, we 
found a column density at the integrated intensity peak position 
of /V(N 2 H + )= 6 + 2 x 10 12 cirT 2 , comparable to other starless 
cores moderately evolved in a chemical sense like e.g. L1512, 
L1498, L1495, L1517B JCrapsi et all l2005h . and a factor of 
2-3 less than in IRAM 041 91, the other very low-luminosity 
Class object known ( Belloche & Andre, 2004). Similarly, the 
N2D + column density derived from the N2D + (l-0) and (2-1) 
IRAM-30m spectra yielded a value of = 6 ± 1 x 10" cm" 2 
from both lines. The deuterium fractionation at the peak of 
L1014 was, thus, found to be 10%, a value higher than the aver- 
age starless core dCrapsi et alll2005l) . N2FL and N 2 D + column 
densities were calculate d also in the Large Velocity Gradient 
approximation (see e.g. ICrapsi et all l2005t for details) yield- 
ing values of 5.3 x 10 12 cirT 2 and 5.2 x 10 11 cirT 2 respectively. 
This approach furnishes also an estimate of the H2 volume 
density from observations of multiple rotational transitions; in 
this way we obtained nfEL;)- 1.9 x 10 5 cirT 3 from N2I-L and 
n(H2)- 4.0 x 10 5 cirT 3 from N2D + , both estimates are consis- 
tent within their errors with the central density evaluated from 
the dust continuum. 

12 CO(1-0) and 13 CO(1-0) profiles are much more asym- 
metric than the other lines and their observed line widths 
(~2 km s _1 and ~1 km s _1 ) are broader (~0.5 km s _1 ). This 
occurrence is not consistent with a line broadening caused 
by opa city effects alone, as already noted bv lRobert & Pa gani 
(1993). A possible explanation could be found considering that 
CO and 13 CO emission comes from the combination of emis- 
sions by randomly-moving, low-density parts of the cloud with 
different velocities along the line of sight and with larger tur- 
bulence. C I7 O(l-0) opacity was found to be 2.4 +1.5 from 
hyperfine structure fitting although the low signal-to-noise ra- 
tio lets us think that this might be just an upper limit. In fact, 
the total column density derived from other indicators is totally 
inconsistent with such a high opacity for C n O. 

Interestingly, our CS(2-1) spectra do not show the asym- 
metric double peaked structure, with the blue peak stronger 
than t h e red one, typical of the c o res undergoing i nfall 
dZhoul Il992l iMardones et ail Il997l Lee. Mvers & Tafallal 
199 9JK This feature was clear ly seen in IRAM 04191 
llAndre. Motte & Bacmannl fl999). Moreover, the difference in 
line velocity between CS and N2H + in units of the N2H + (1- 
0) line width (5Vcs = (Vcs - Vn 2 h+)/AVn 2 h+X that is sup- 
posed to gauge the extended infall motions of the core 



dMardones et all Il997t iLee. Mvers & Tafallal Il999tl . is very 
small inL1014(0.04). 

On the contrary, HCO + (3-2) observations show the oppo- 
site asymmetry revealing a red peak brighter than the blue peak 
(see Figure|3}. This line profile can be produced by a rotating, 
expanding, or pulsating core jLada et all 12003). We note that 
the velocity of the self absorbed feature corresponds to the ve- 
locity of the optically thin C I8 0(2-1) line as predicted by the 
models. Given the moderate signal to noise ratio of the present 
data, it is necessary to confirm these indications with higher 
sensitivity and higher resolution observations and to retrieve 
a map of HCO + in order to search for spatial patterns for the 
expansion asymmetry. 

We failed to detect H 13 CO + (3-2) and HCN(3-2) towards 
L1014-IRS up to a 0.02 and 0.03 K km s" 1 level respectively. 

3.3. The search for the outflow in L1014 

In Figure [6j we present channel maps of CO(1-0) (top) 
and 13 CO(1-0) (bottom) in a range of velocities, -9 km s _1 
to 15 km s _1 , centered around the L1014 typical velocity 
(~4.2 km s _1 ). Emission was integrated in 2 km s _1 intervals. 
No CO emission was found anywhere in the mapped area be- 
tween -35 km s _1 and -9 km s _1 and between 15 km s _1 and 
25 km s _1 , i.e at velocities departing more than 10 km s _1 from 
the L1014 rest velocity. The Icr errors for the 12 CO and 13 CO 
channel maps are 0.06 K km s _1 and 0.03 K km s _1 respec- 
tively, less than the lowest contours adopted in Figures [6] and 

m 

Using IRAM 04191 dAndre. Mott e & Bac manrlfl999T) as a 
guide, we searched for 12 CO(1-0) wing emission at velocities 
3 km s _1 < |V — Vol < 9 km s -1 . No wing-like emission was 
found in this range (see maps labeled -5, -3, -1 and 7, 9, 11 
in Figure|6j. 

Nevertheless, given the very low luminosity of L1014-IRS 
we examined also the possibility of less energetic wings. The 
most likely outflow feature present in L1014 is along position 
angle 45° (measured East of North) where we can see red emis- 
sion in the SW lobe between 5 and 7 km s and blue emission 
in the NE part in the bin 1 to 3 km s -1 . This axis was identified 
mainly from the jet-like appearance of the "red lobe", but the 
fact that this feature is also seen in the 13 CO(1-0) and C ls O(l- 

0) emission, which typically are not sensitive to outflows, poses 
serious questions as to whether this is indeed due to a classical 
outflow lobe. Moreover, as can be seen in the sequence of spec- 
tra along this axis (Figure^, we do not see the typical outflow 
features in the line profiles. For example: 

1) 12 CO(1-0) line at the cloud rest velocity is not present 

throughout all the cut and is not the dominant component 
in all the positions, 

ii) the "low velocity wings", seen for example in the positions 
(-500, -500) in the red lobe or (100,100) in the blue lobe, 
look less like wings in the following positions (-200, -200) 
and (400, 400). 

From the spectra in Figure we deduce that the red lobe and 
the blue extended emission seen in the 12 CO map arise from 
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different parts of the cloud with different velocities along the 
line of sight. This view is reinforced also by the observations of 
13 CO and C ls O emission in the same region and with consistent 
velocities. The presence of gas with different velocities along 
the same line of sight would also explain the relatively broad 
12 CO and 13 CO line w idths found towards the entire map (see 
iRobert & Paganil 19931 and Sectionl3~2l 

A peculiar feature of the 12 CO channel maps is the pres- 
ence of several small spots of emission in the NW part of the 
map at velocities smaller than 1 km s _1 and one spot at 9 km 
s _1 in the SE. Although these may be interpreted as molecu- 
lar "bullets" coming from the L1014 nucleus, their asymmetric 
positioning, low velocity and large distance from the nucleus 
(0.0 5 pc) make th em very different from the prototypical cases 
iBachmerlll996l) . 

We conclude that no classical outflow signature is present at 
the large scales investigated with the FCRAO. Further sensitive 
observations must be performed with smaller beams to probe 
the CO emission within the inner 45" (the FCRAO beam size). 

3.4. The background component 

The Spitzer observations presented bv lYoung et all fc004 left 
open the possibility of L1014-IRS being a more massive but 
less embedded distant protostar aligned by chance with the 
L1014 nucleus. Although the relative rarity of embedded proto- 
star and starless cores do not favour this chance, the fact that we 
are looking close to the galactic plane and towards the Perseus 
spiral arm did not allow us to disregard it just on a statistical 
basis. Now, recent deep near-infrared observations confirm that 
L1014-IRS is associated with the nearby L1014 core (T. Huard 
et al. 2005, in prep.), our FCRAO observations are consistent 
with this finding. 

We searched for signatures of infall-outfiow activity to- 
wards L1014-IRS also at the Perseus arm velocity (near - 
40 km s- 1 ). No N 2 H + (l-0), CS(2-1) and C 18 O(l-0) emis- 
sion was detected at these velocities, most likely because of 
beam dilution given the small scales traced b y those species 
and t h e large Perseus arm distance (2.6 kpc. iBrand & Blitzl 
ll993l) . lYoung et al.N2004 calculated that, assuming a distance 
of 2.6 kpc for L1014-IRS, the protostar whose emission would 
fit best the observed SED would have a luminosity of 16 L Q . 
Class I and Class protostars with com parable bolometric lu- 
minosity li ke LI 165 (TVisser Richer & Chandlerll2002l) . L1448 
and LI 157 jBachillerll 19961) drive outflows of ^0.3 pc in size, 
corresponding to ^20" at the Perseus arm distance. Although 
this size is smaller than the FCRAO beam, we performed a 
search for high velocity wings in the -40 km s _1 component. 
CO and 13 CO channel maps at these velocities are shown in 
Figure^ with the same intensity scale as in Figure^ No emis- 
sion was detected below -58 km s and above -36 km s . 
As for the component around 4 km s _1 , no clear-cut outflow 
features are seen in the channel map. 

4. Discussion 

The contradiction between the relatively young evolutionary 
status of L1014 derived from chemical indicators and the pres- 



ence of a protostellar object embedded in its nucleus is puz- 
zling. L1014 shows continuum and line intensities much lower 
than well-studied starless cores; in particular, the H2 volume 
density, N2H + column density and the degree of CO depletion 
are smaller than the average starless core, from which we did 
not expect L1014 to be close to the star formation. Moreover, 
the CS observation does not indicate the presence of inward 
motions, and the HCO + line profile show hints of outward mo- 
tions. 

On the other hand, we report also relatively high deuterium 
fractionation and broad (if compared to other low-mass starless 
cores) N2H + and N 2 D + lines which are typical of more evolved 
objects. 

We speculate that either i) the v ery low-luminosity 
nature of the central source (0.09 L Q , Young etaU l2004 
makes its chemical and dynamical evolution different 
to any other observed core or ii) that a low-mass, low- 
luminosity "seed" is present lo ng before the protos tellar 
phase in every "starless" cor e. iBoss & Yorka \ 19951) and 
Mas unaga. Mivama & Inutsukal Jl998l) modeled the early 
phase of protostar formation introducing the "first hydrostatic 
core" or "Class -I" protostar as a short lived precursor of the 
Class 0. During this phase, the central temperature should 
reach a val ue of 200 K (3 ti mes lower than the temperature 
derived by lYoung et alJ 120041) and no outflow emission is 
expected. According to the duration of this phase, we could 
expect the detection of similar sources in the nuclei of some 
other starless cores in future Spitzer observations. 

Consid ering that L1 14-IR S was classified as a Class 
protostar jYoung et all 120041) having met the require- 
ments of Thni < 70 K and Lh n \J L smm < 200 
llAndre. Ward-Thompson. & Barsonvi Il993h . our observations 
indicate that L1014-IRS might be the first Class proto- 
star not associated with an observed molecular outflow. Even 
IRAM 04191, which has a bolometric luminosity comparable 
to L1014-IRS (~0. 15 Lp), has an easily det ected, very extended 
outflow (And re. Motte & BacmannL Il999l) : thus, it seems un- 
likely that the weakness of the central source is responsible 
for the lack of the outflow detection. In any case, one should 
bear in mind that differences in the inclination angle, collima- 
tion factor, in the age, as well as in the external environment 
could make the detection of the outflow more difficult even in 
presence of two jets with comparable momenta. We can spec- 
ulate that either the accretion rate of L1014-IRS is currently 
too little to power the outflo w or that the mag netic field that 
threads the disk is too weak (Bachiiierl fl99fl) . Alternatively, 
we may have caught L1014 in the epoch betw een the proto- 
star formation and the outflow ignition (see e.g. lBoss & Yorkd 
1995; Masunaga. Mivama & Inutsuka 1998), although this hy- 
pothesis seems less likely given the short lifetime associated to 
this phase. 

We remark that our observations do not probe small scale 
outflows that would be diluted in the 45" FCRAO beam, in 
particular if the outflow is very young and therefore compact. 
Strong molecu lar outflows are a hallmark of Class protostar 
llAndrel fl994) , thus it is of fundamental importance to extend 
the search for the outflow to a smaller scale. 
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The present data alone cannot rule out conclusively the 
possibility of chance alignment, although our search for back- 
ground dense cores through N2H + and CS was negative, and no 
signatures of molecular outflows were found in the CO back- 
ground component. The absence of a background dense core 
could be explained by the chance crossing of a T-Tauri star 
born elsewhere and expell ed in the direction o f L1014 (see 
the case of PV Cephei, Goodm an & Arcell2004h . This expla- 
nation, although difficult to rule out, is even less probable than 
the chance alignment, given the additional requirements of a 
favourable trajectory from the parent cloud to L1014 and the 
coincidence in the epoch of observations. 
Further evidence that L1014-IRS is a young stellar object em- 
bedded within L1014 comes from deep near-infrared observa- 
tions and will be presented by T. Huard et al. (in preparation). 

5. Conclusions 

We observed the starless core L1014 with the FCRAO an- 
tenna in N 2 H + (l-0), CS(2-1), C 18 O(l-0), 13 CO(1-0) and 
12 CO(1-0), combined with literature 1.2-mm continuum data 
and new N 2 H + (l-0), N 2 H + (3-2), N 2 D + (l-0), N 2 D + (2-l), 
N 2 D + (3-2) and C 17 O(l-0) observations from the IRAM-30m 
and C 18 0(2-l), C n O(2-l), HCO + (3-2), H 13 CO + (3-2) and 
HCN(3-2) spectra from CSO to study its chemical status and 
to search for the presence of a molecular outflow. The results 
of our study are summarized below. 

1. The chemical and physical properties of L1014 derived 
from the present observations are not typical of highly evolved 
low-mass starless cores. In particular, we found: molecular hy- 
drogen volume density of n(H 2 )^ 2.5xl0 5 cirT 3 , N 2 H + column 
density of JV(N 2 H + )^ 6 x 10 12 cirT 2 , M vir = 2.1 M , CO inte- 
grated depletion factor equal to 9, absence of CS double peaked 
profile with infall signature, absence of velocity shifts between 
self-absorbed optically thick (CS) and optically thin (N 2 H + ) 
tracers and a profile asymmetry in HCO + (3-2) consistent with 
outward motions. On the other hand, we found an enhanced de- 
gree of deuterium fractionation equal to 10% which is higher 
than the average starless core, and broad N 2 H + and N 2 D + lines 
suggestive of unresolved kinematical activity in the inner nu- 
cleus. 

These diverging indications could be reconciled consider- 
ing that we are observing either a very young stage of star 
formation or an extremely low-luminosity object exhibiting 
characteristics that differ significantly from previously known 
cases. 

2. No classical signatures of a molecular outflow were 
found towards L1014. In particular, we note the absence of high 
velocity wing and of symmetric well-defined red-blue lobes in 
the CO channel maps. Bearing in mind that the scale probed 
by the present observations might be too large to detect the 
outflow in L 1014, we suggest that the formation of a protostar 
might occur prior to or in absence of the molecular outflow. 

The presence of a protostar in L1014 seems to challenge 
the idea that there is a unique path to forming a protostar that 
all cores must follow. 
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Table 1. FCRAO settings and parameters. 



line 


frequency 


HPBW 




Av res 


vel. cov. 


Mapped region 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


N 2 H + (l-0) 


93.1737725 


54 


280 


0.157 


-94 -> 57 


5'x 5' 


CS(2-1) 


97.980953 


52 


280 


0.149 


-90 -» 55 


5'x5' 


C 18 O(l-0) 


109.782173 


46 


280 


0.133 


-60 -» 70 


5'x 5' 


13 CO(1-0) 


110.201370 


46 


440 


0.133 


-81 -> 45 


20' x 32' 


12 CO(1-0) 


115.271203 


44 


920 


0.127 


-78 -» 43 


20' x 32' 



Note. - Col. (2) line rest frequency (GHz); Col. (3) Half Power Beam Width ("); 

Col. (4) System Temperature (K, main beam scale); Col. (5) Channel Spacing (km s -1 ); 

Col. (6) Usable velocity coverage (km s -1 ); Col. (7) R.A. and Dec. extension of the mapped region. 



Table 2. IRAM-30m and CSO telescope settings and parameters. 



line 


frequency 


HPBW 




Av res 


vel. cov. 


N obs 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


IRAM-30m observations 


N 2 H + (l-0) 


93.1737725 


26 


165 


0.021 


-10 -> 16 


5 


N 2 H + (3-2) 


279.511863 


9 


1500 


0.042 


-16 -» 12 


3 


N 2 D + (l-0) 


77.109626 


32 


155 


0.025 


-14 -» 20 


5 


N 2 D + (2-l) 


154.217137 


16 


300 


0.026 


-10 -> 11 


3 


N 2 D + (3-2) 


231.321966 


11 


440 


0.050 


-16 -> 16 


5 


C 17 O(l-0) 


112.358990 


22 


240 


0.034 


-10 -> 16 


1 


CSO observations 


C 18 0(2-l) 


219.560352 


34 


390 


0.130 


-29 -» 37 


13 


C 17 0(2-l) 


224.714214 


33 


330 


0.128 


-29 -» 37 


5 


HCO + (3-2) 


267.557620 


28 


320 


0.107 


-21 -» 28 


1 


H 13 CO + (3-2) 


260.255478 


29 


300 


0.110 


-24 -> 32 


1 


HCN (3-2) 


265.886434 


28 


500 


0.054 


-24 -» 31 


1 



Note. - Col. (2) line rest frequency (GHz); Col. (3) Half Power Beam Width (") 

Col. (4) System Temperature (K, main beam scale); Col. (5) Channel Spacing (km s -1 ); 

Col. (6) Usable velocity coverage (km s -1 ); Col. (7) Number of observed positions. 
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Table 3. Line parameters at the 1 .2-mm continuum peak position from line profile fitting. 
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line 


intensity 


Vlsr 


AV 


T 


(1) 


(2) 


(3) 


(4) 


(5) 


FCRAO spectra 


N 2 H + (l-0) 


0.52+0.03 


4.239+0.014 


0.354+0.033 


4.6+3.2 


CS(2-1) 


0.41+0.02 


4.279+0.014 


0.677+0.031 




C 18 O(l-0) 


0.84+0.02 


4.224+0.005 


0.484+0.012 




13 CO(1-0) 


3.18±0.07 


4.258+0.008 


0.864+0.019 




12 CO(1-0) 


9.30+0.17 


4.211+0.022 


2.261+0.049 




IRAM-30m spectra 


N 2 H + (l-0) 


1.850+0.032 


4.242+0.004 


0.354+0.009 


6.1+0.8 


N 2 H + (3-2) 


0.263+0.036 


4.270+0.035 


0.443+0.058 


<0.100 


N 2 D + (l-0) 


0.203+0.010 


4.248+0.007 


0.284+0.020 


2.1+1.6 


N 2 D + (2-l) 


0.182+0.015 


4.273+0.010 


0.307+0.029 


<0.100 


N 2 D+(3-2) 


0.066+0.014 


4.303+0.044 


0.315+0.128 


<0.100 


C 17 O(l-0) 


0.476+0.022 


4.229+0.009 


0.310+0.027 


2.4+1.5 


CSO spectra 


C 18 0(2-l) 


0.641+0.036 


4.270+0.019 


0.575+0.048 


1.4° 


C 17 0(2-l) 


0.344+0.035 


4.092+0.023 


0.489+0.051 


< 0.100 


HCO + (3-2) 


0.316+0.015 


4.410+0.019 


0.819+0.041 




H 13 CO + (3-2) 


0.025+0.019 








HCN (3-2) c 


<0.03 









we integrated over all the components; 12 CO and 13 CO lines were 

integrated between 1 and 7 km s -1 and 3 and 5.5 km s -1 respectively); 

Col. (3) Rest velocity (km s _I ) 

Col. (4) Full width half maximum (km s ); 

Col. (5) Sum of the opacity of all the hyperfine components. 

a Opacity was evaluated from the opacity of C 17 (2-1) and assuming the 

relative abundance as in Wilson & Rood (1994). 

c Upper limit on integrated intensity was evaluated for a line-width 

equal to H 13 CO+(3-2). 



Fig. 1. FIGURE 1 IS ATTACHED AS A JPEG FILE IN ASTRO-PF ARCHIVE. Optical image of the 27'x37' field around L1014 
from the Digita l Sky Survey. Overlaid onto the ima ge are contours of beam-averaged visual extinction, following the NICE 
technique (e.g.. lLada et all Il994l lAlvesetall F998) of convolving line-of-sight measurements of the H-K color excesses of 
2MASS sources with a 30" Gaussian beam. The contours are drawn at Ay=[3, 5, 7, 9, 11] magnitudes. Our FCRAO observations 
mapped much of this region in CO(l-O) and 13 CO(1-0), while the N 2 H + (l-0), CS(2-1) and C 18 O(l-0) observations covered 
Mb area represented by the rectangle in the €gW&. <Eii01*l4§ls%giPl?P \fer efeflfteYaSfJffie field, while B362 can be seen 10' north of 
L1014. The cross indicates the position of L1014-IRS, the point source observed with the Spitzer Space Telescope consistent 
with being an embedded protostar. 




100 -100 -200 -300 100 -100 -200 -300 

Aa (") Aa (") 

Fig. 2. 1.2-mm continuum, N2H + (l-0), CS(2-1) and C I8 O(l-0) emission towards L1014. Molecular line maps were taken at 
FCRAO with the exception of the N 2 H + (l-0) map in the inset (fro m IRAM-30m) and the C ls O(2-l) map in the inset (from 
CSO) The 1.2-mm continuum was taken from lYoung et alJ J2004) and it was observed with IRAM-30m. Contour levels start 
at and increase by 0.08 K km s for the FCRAO molecular data and 4 mJy/(ll" beam) -1 for the dust continuum . The white 
contour represents the half peak intensity in all the maps. Beam sizes are displayed at the top-left of each map but the 1.2-mm 
map. The white cross places the position of a point like source detected by the Spitzer Space Telescope with colors compatible 
with an embedded protostar. Angular offsets are relative to the 1.2-mm continuum peak (21 h 24 m 07:5, 49°59'05", J2000). The 
inset on the N 2 H + (l-0) FCRAO map shows the 5 points observations of N 2 H + (l-0) at IRAM-30m. Levels start at 0.8 K km 
s~'and increase by 0.3 K km s . The IRAM-30m beam size is shown inside that of the FCRAO. Similarly, the inset in the 
C 18 O(l-0) FCRAO map shows the 13 points map of C ls O(2-l) taken at CSO. Levels stall at 0.35 K km s 'and increase by 0.1 K 
km s . The CSO beam size is shown inside that of the FCRAO. 
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Velocity (km s ) 

Fig.3. CO(l-O), I3 CO(1-0), C 18 O(l-0), CS(2-1) and N 2 H + (l-0) spectra observed with FCRAO towards the 1.2-mm peak of 
L1014. Gaussian (or hyperfine in the case of N2H + (l-0)) fits are plotted. 
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Fig. 4. N 2 H + (l-0), N 2 H + (3-2), C I7 O(l-0), N 2 D + (l-0), N 2 D + (2-l) and N 2 D + (3-2) spectra observed with the IRAM-30m to- 
wards the 1.2-mm peak of L1014. Hyperfine fits are plotted. 
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H 13 CO + (3-2) 




2 4 6 8 

Velocity (km s ) 

Fig. 5. C 18 0(2-l), C n O(2-l), HCO + (3-2), H 13 CO + (3-2) and HCN(3-2) spectra observed with the CSO towards the 1.2-mm 
peak of L1014. Gaussian or hyperfine fits are plotted. The asymmetric shape of HCO + (3-2) is consistent with outward motions 
of the gas. 
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Fig. 6. Channel maps of the CO(l-O) (top panel) and 13 CO(l-0)(bottom panel) lines around the rest velocity of L1014. Contours 
levels start and increase by 0.8 K km s _1 for the CO data and by 0.05 K km s for the 13 CO ones. Velocity bins are spaced 
by 2 km s'and the starting velocity of the bin is reported in the top right of each map. The central cross in each map indicates 
the dust peak while the other four delineates the area shown in Figure|5] The solid line shows the direction of the cut studied in 
Figure 
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Upper bold line: C0(1— 0) 
Lower light line: 13 CO(1-0) (X 3) 




5 10 5 10 

Velocity (km/s) Velocity (km/s) 



Fig. 7. CO(l-O) (upper bold line) and 13 CO(1-0) (lower light line) spectra along the direction of the best-candidate outflow. 
South-West (red) branch is in the left panel and North-East (blue) branch is in the right panel. I3 CO(1-0) spectra were multiplied 
by 3. 
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Fig. 8. Channel maps of the CO(1-0) (top panel) and 13 CO(l-0)(bottom panel) lines for a second feature visible at ~ -40km s 1 
and probably associated with the Perseus Arm. Contours levels start and increase by 0.8 K km s 1 for the CO data and by 0.05 K 
km s _1 for the I3 CO ones. Velocity bins are spaced by 2 km s~'and the starting velocity of the bin is reported in the top right of 
each map. The central cross in each map indicates the dust peak while the other four delineates the area shown in Figure|5] 
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